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Analyses are carried out to assess the impact of nuclear data uncertainties on k e f / for the European 
Lead Cooled Training Reactor (ELECTRA) using the Total Monte Carlo method. A large number of 
Pu-239 random ENDF-formated libraries generated using the TALYS based system were processed 
into ACE format with NJOY99.336 code and used as input into the Serpent Monte Carlo neutron 
transport code to obtain distribution in k e ff. The fc e // distribution obtained was compared with 
the latest major nuclear data libraries JEFF-3.1.2, ENDF/B-VII.l and JENDL-4.0. A method 
is proposed for the selection of benchmarks for specific applications using the Total Monte Carlo 
approach. Finally, an accept/reject criterion was investigated based on \ 2 values obtained using 
the Pu-239 Jezebel criticality benchmark. It was observed that nuclear data uncertainties in fe e // 
were reduced considerably from 748 to 443 pcm by applying a more rigid acceptance criteria for 
accepting random files. 



I. INTRODUCTION 

In recent times, there has been an increased use of nu- 
clear reaction codes in the nuclear data evaluation pro- 
cess for the computation of cross sections, spectra and 
angular distributions required for a large variety of ap- 
plications pQ. The use of model codes in the evaluation 
of nuclear data offers several advantages, among them; 
the preservation of energy balance and coherence of par- 
tial cross sections together with total and reaction cross 
sections, the prediction of data for unstable nuclei and 
providing data where experimental data are unavailable. 
Where experimental data are available, they are used for 
fine tuning input parameters to model codes. The Total 
Monte Carlo (TMC) method, developed around TALYS 
by the Nuclear Research Group, Petten, has the capa- 
bility of incorporating microscopic nuclear physics and 
macroscopic nuclear reactor design into one simulation 
scheme [5] . In this way, uncertainties can be assigned to 
macroscopic reactor parameters together with their sen- 
sitivity to uncertainties in specific cross-sections. The 
probability distribution for various reactor parameters 
can also be observed. TMC has been tested for a va- 
riety of systems [3]. In this paper, we applied the Total 
Monte Carlo approach in assessing the impact of Pu-239 
cross section uncertainties on the full core 3-D SERPENT 
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Monte Carlo model of the European Lead Cooled Train- 
ing reactor(ELECTRA). ELECTRA is a 0.5 MW reactor 
proposed within the GENIUS (Generation IV research in 
universities of Sweden) framework [3] . It aims at devel- 
oping technology necessary for safe and economic deploy- 
ment of Generation IV nuclear reactors, in particular the 
lead cooled fast reactor. Since Pu-239 makes up of 30% of 
the total mass of the fuel, its uncertainties are expected 
to impact significantly on the core behavior of the reac- 
tor. A novel approach for the selection of benchmarks 
for applications is present in this work. The selection of 
benchmarks for applications is usually done with 'by eye' 
judgment. This method is dependent on the experience 
of the user and therefore not suitable for the Total Monte 
Carlo methodology which lays emphasis on automation 
and reproducibility [2]. Also, an accept/reject criterion 
for reducing uncertainties in reactor macroscopic param- 
eters based on limiting % 2 values is investigated using the 
Pu-239 Jezebel criticality benchmark. 



II. METHODOLOGY 

A. Simulations 

The input files used in this study are the SERPENT 
geometry input file developed at KTH, Sweden [J] and 
740 random ENDF files obtained from the TENDL 
project 0. The Pu-239 Jezebel benchmark MCNP in- 
put file obtained from the International Handbook of 
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Evaluated Criticality Safety Benchmark Experiments [5] 
was converted into a SERPENT file for simulations. The 
random files were processed into ACE format using the 
Njoy99.336 processing code at 1200K. Simulations were 
performed with the Serpent Monte Carlo code version 
1.1.17 6J for the ELECTRA core at zero burnup with 
the absorber drums set at startup position. Criticality 
calculations were performed for 20 inactive and 500 k e f / 
active cycles with 50,000 neutrons per cycle correspond- 
ing to 25 million histories resulting in a statistical uncer- 
tainty in k e ff of 21.5 pcm. ENDF-6 formatted files for 
Pu-239 obtained from the latest major nuclear data li- 
braries JEFF-3.1.2, ENDF/B-VII.l and JENDL-4.0 were 
processed in a similar way and criticality calculations per- 
formed. The results obtained are compared with mean of 
the k e ff distribution of the random nuclear data files. 

B. Fast TMC 

The fast TMC methology presented in [7] is used in 
this work. However, instead of using it to obtain short 
computation times it is used to obtain a very accurate and 
precise estimate of the uncertainty in reactor parameters 
due to nuclear data ((t nd ). In TMC, (<J ND ) is obtained 
by running N different ND files. The result is a spread in 
the data (<7 & s ) that is both due to statistics (<J s tat) and 
due to nuclear data (cr ND ): 



Since both the parameter, e.g, fc e //j and its statistical un- 
certainty, Ak e ff, is determined by most transport codes, 
we obtain a distribution of both k e f t and Ak e f f when we 
run with our N random files. The spread in the distribu- 
tion of k e f f is then <j bs and the Ak e f / is equal to a s t a t ■ 
Eq. 1 holds assuming that there is no correlation between 
statistics and nuclear data; the correlation between k e ff 
and the corresponding statistical error was investigated 
and no correlation was observed. The methodology pro- 
duces an accurate estimate of a ND if Ak e ff is accurately 
determined by the transport code. To test this we exe- 
cuted the serpent input file 740 times with constant nu- 
clear data, but with different random seeds. If a b s = 
Ak e ff then the transport code determines Ak e ff cor- 
rectly; in our case with SERPENT and ELECTRA and 
this was found to be true. 

C. Selecting benchmarks for applications 

Criticality benchmarks are used mostly for validation 
of calculational techniques in reactor calculations and for 
establishing minimum subcritical margins for operations 
with fissile material [5], In the nuclear data evaluation 
process, benchmarking evaluated data against differen- 
tial and integral data is very important. This is mostly 
done by selecting large sets of criticality, dosimetry, fu- 
sion, activation benchmarks among others. The selec- 
tion varies from one application group to another as a 



results of different expertise, purpose of the evaluation 
and accessibility to benchmarks [3]. To solve this prob- 
lem, we proposed the TMC methodology for the selec- 
tion of benchmarks. First, random files obtained from 
the TENDL project are processed into usable format for 
reactor codes. Criticality simulations are performed for 
1) the i th benchmark under consideration and for 2) the 
specific reactor system using same processed random files. 
A correlation test is then performed by plotting the fc e / / 
obtained for the reactor system against that of the bench- 
mark. If a strong correlation exists between the two sys- 
tems, this can be interpreted as a good representation of 
your system. This methodology has been tested only on 
the Pu-239 Jezebel benchmark and we plan to investigate 
this in more detail taking into consideration fusion and 
shielding benchmarks. The proposed method also opens 
up possibilities for assigning weights to each benchmark 
based on the goodness of fit or the correlation coefficient. 
In Fig. 1, a correlation between the k e ff computed for 
ELECTRA and that of the Pu-239 Jezebel benchmark is 
presented. 
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FIG. 1. Correlation between fc e // determined for ELECTRA 
and the Pu-239 Jezebel benchmark for Pu-239 random files. 



D. Accept /reject criteria based on benchmarks 

An accept /reject criterion is investigated to study how 
limiting \ 2 values for accepting random files could reduce 
nuclear data uncertainties in safety parameters. Simula- 
tions were performed for Pu-239 random files [2] and 
with Pu-239 Jezebel criticality benchmark and with the 
ELECTRA Serpent Monte Carlo code input file. Limit- 
ing chi squares were set and corresponding nuclear data 
uncertainty calculated. In Fig. 2, we present a plot of lim- 
iting chi square verses Pu-239 random files index. 

III. RESULTS 

A distribution in k e ff was obtained for 740 Pu-239 
random ENDF files and the mean compared with results 
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from other nuclear data libraries. The JEFF-3.1 (the ref- 
erence library) temperature dependent cross-section li- 
braries were maintained for all other isotopes except for 
Pu-239 files which was changed with files from TENDL, 
JENDL-3.1.2, JENDL-4.0 and END/B-VII.l. It was 
observed that the mean value of of the random files 
1.001540 ± 0.00021 compared favorably with the other 
nuclear data libraries: ENDF/B-VII.l was 1.00067 ± 
0.00021, JEFF-3.1. 1 was 1.00023±00022, JENDL-4.0 was 
1.00043 ± 0021. This can be attributed to the contri- 
bution of the well known fission cross section of Pu-239 
which has a relatively high impact on the effective multi- 
plication factor. The uncertainty due to nuclear data was 
determined and found to be 745 ± 19 pcm. In Fig. 2, we 
present a plot of limiting x 2 values against the random 
file index. Nuclear data uncertainties were determined 
for each limit and the results is illustrated in Fig. 3 
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FIG. 2. Limiting chi square against Pu-239 random file index. 

The correlation between k e f / computed for ELECTRA 
and for the Pu-239 Jezebel benchmark is shown in Fig. 1. 
It can be seen from the plot that there exists a strong 
correlation between the two systems implying that, the 
Pu-239 Jezebel benchmark is a good representation of 
ELECTRA. The Pearson correlation coefficient which is 
a measure of the strength of the relationship between the 
two systems was computed and found to be R=0.84, fur- 
ther underscoring the strong correlation observed. For 
the computation of \ 2 f°r a large set of benchmarks, 
we propose that a weighted factor (uii), which should 
be equal to the absolute value of the correlation coeffi- 



cient (Ri) be set for the i th benchmark. We plan to test 
this methodology on a large set of benchmarks and for 
different nuclides. The approach has however been ap- 
plied Am-241 in burnup evaluation of ELECTRA [3]. 
It was observed that, a 25% reduction in nuclear data 
uncertainty was achieved after a setting a more rigid ac- 
ceptance criteria. In Fig. 3, we present the nuclear data 
uncertainty against limiting x 2 • it can be observed from 
the figure that by setting a minimum limiting x 2 of 5.31e- 
06, we were able to reduce the uncertainty in nuclear data 
from 748 to 443 pcm. This therefore implies that, by im- 
plementing an accept /reject criterion in the Total Monte 
Carlo method for accepting random files, nuclear data 
uncertainty can be reduced significantly. 




Limiting chi squares x io~ 4 

FIG. 3. Nuclear data uncertainty against limiting chi square 
using the Pu-239 Jezebel benchmark and Pu-239 random files. 



IV. CONCLUSIONS 

A study has been carried out by using benchmarks 
within the Total Monte Carlo methodology on nuclear 
data uncertainty propagation for the European Lead 
Cooled Training Reactor (ELECTRA). It was observed 
that, nuclear data uncertainties could be reduced consid- 
erably in criticality calculations after introducing an ac- 
cept/reject criterion based on integral benchmarks into 
the TMC calculation chain. A method for the selection 
of benchmarks with TMC was also proposed which opens 
up possibilities for assigning weights to benchmarks. 
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